Cross section, final spin and zoom-whirl behavior in high-energy black hole collisions 
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We study the collision of two highly boosted equal-mass, nonrotating black holes with generic 
impact parameter. We find such systems to exhibit zoom-whirl behavior when fine tuning the 
impact parameter. Near the threshold of immediate merger the remnant black hole Kerr parameter 
can be near maximal (o/M > 0.95) and the radiated energy can be as large as 35 ± 5% of the 
center-of-mass energy. 
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J. Introduction. As the major foundation of experi- 
mental particle physics, high energy collisions enable us 
to probe the fundamental characteristics of short range 
interactions: the larger the center-of-mass (CM) energy 
of the collision, the shorter the distance one can probe. 
The mass-energy equivalence of relativity, however, im- 
poses a distinct limit on this method. All forms of en- 
ergy gravitate and collisions at trans-Planckian energies 
are expected to generically form black holes (BHs). At- 
tempts at probing shorter distances by increasing the CM 
energy of the system are in consequence barred by the 
formation of the event horizon Because the gravi- 
tational field in this regime is predominantly sourced by 
the kinetic energy, the nature of the interaction should be 
rather insensitive to the internal structure of the parti- 
cles Q . The trans-Planckian scattering of point particles 
should therefore be well described by BH scattering Q . 

Ultra-relativistic BH scattering simulations are rele- 
vant for recent proposals to solve the hierarchy problem 
by adding "large" extra dimensions [3] or an extra dimen- 
sion with a warp factor ||, which allow for the Planck 
scale to be near the electroweak scale. This offers the 
exciting possibility to produce BHs in particle colliders 
and ultra high-energy cosmic ray interactions with the at- 
mosphere 0, ||| . Classical BH production cross sections 
and the energy radiated in gravitational waves (GWs) 
are inputs for the Monte-Carlo simulators used to search 
for such events @- Further applications of high-speed 
BH collisions to high-energy physics have recently been 
suggested by the AdS/CFT correspondence [1, @. 

This Letter reports on several remarkable features aris- 
ing in ultra-relativistic BH collisions with generic impact 
parameter, thus generalizing the head-on collision results 
of We also extend a recent preliminary exploration 
by Shibata et al. [ll|, which estimated the impact pa- 
rameter at the threshold between BH formation and scat- 
tering for different boost magnitudes. Unless stated oth- 
erwise, we use geometrical units G — c = 1. 
II. Setup. The simulations presented in this work have 
been obtained with the Lean code il2l. which is based 
on the Cactus computational toolkit [13| . The code em- 
ploys mesh refinement (provided by Carpet [3]) and 
the apparent horizon (AH) finder AHFinderDirect 



jl5| . Puncture initial data are provided by a spectral 
solver 11611. The simulations are performed as in the head- 
on case [10| , but here we use equatorial instead of octant 
symmetry and employ higher resolution near the BHs. 

We set up a coordinate system such that the BHs start 
on the x-axis separated by a coordinate distance d and 
with radial (tangential) momentum P x {P y ). The impact 
parameter is b = L/P = P y d/P, where P is the linear mo- 
mentum of either BH, and L is the initial orbital angular 
momentum. We extract gravitational radiation by com- 
puting the Newman-Penrose scalar '5 4 at different radii 
r ex from the center of the collision. '5 4 is decomposed into 
multipoles %/)i m using spin- weight —2 spherical harmonics: 



1=2 2^m=- 



iYi m (0 ,(/>) ipi m (,t,r c 



where 9 is measured relative to the z-axis. The esti- 
mated spurious or "junk" radiation in the initial data is 
quite insensitive to the impact parameter and compara- 
ble to that present in the head-on case ; we remove it 
from reported results in the same manner. Errors due to 
discretization and finite extraction radius are compara- 
ble to those reported in [l(| • The estimated uncertainties 
in radiated quantities are 3% and 15% for low and high 
boost, respectively, and the phase error in the GW signal 
used in the zoom-whirl analysis below is 0.2 rad. 

Our analysis of grazing collisions is based on three one- 
parameter sequences of numerical simulations of equal- 
mass, nonspinning BH binaries. Each sequence is char- 
acterized by fixed initial coordinate separation d and 
Lorentz boost 7 = (1 — u 2 ) -1 / 2 of the holes, while the 
impact parameter b is varied. These sequences arc (1) 
7 = 1.520 (v = 0.753) and d/M = 174.1: (2) 7 = 1.520 
(v = 0.753) and d/M = 62.4; (3) 7 = 2.933 (v = 0.940) 
and d/M = 23.1, where M is the total BH mass. 

III. Results. The results of our study are most conve- 
niently presented in terms of three distinct regimes we 
encounter as the impact parameter is increased starting 
from the head-on limit 6 = 0: (i) immediate mergers, (ii) 
nonprompt mergers and (iii) the scattering regime where 
no common AH forms. These regimes are separated by 
two special values of b: the threshold of immediate merger 
b* and the scattering threshold 6 sca t. The remarkable fea- 
tures of BH binaries in these different regimes will be 
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described in detail in the remainder of this section. 

The scattering threshold 6 scat is defined such that the 
two BHs merge for b < 6 scat and scatter to infinity for b > 
frscat- By analyzing collisions with CM velocity v < 0.90, 



Shibata et al. estimate 6 scat /M ~ (2.5 ± 0.05) /v. 



The analysis of sequences 1 and 2 shows that merger oc- 
curs only for b sca x/M < 3.4, consistent with [ll|, but for 
sequence 3 we find 2.3 < b sca t/M < 2.4, indicating that 
[111 ] may overestimate 6 sca t for large 7. Previous studies 
in the literature, which are based on the Penrose con- 
struction, look for AHs in the union of two shock waves 
and find b scat /M = 1.685M as v -» 1 [l7j. Our results 
suggest that estimates of BH production cross sections 
(oc 6g Cat ) obtained through that construction are accu- 
rate to within a factor < 2. We further emphasize the 
surprising agreement between our simulations and the 
point-particle approximation: for example, a cross sec- 
tion estimate from high-energy scattering off a Kerr BH 
with j ~ 0.98 gives b scat /M ~ 2.36 [jj]. 

The above definition of the scattering threshold is 
purely based on the nature of the end state of the bi- 
nary but ignores details of the interaction. A closer look 
at these details reveals the existence of a threshold of im- 
mediate merger b* [l9l ]. Roughly speaking, for b < b* 
merger occurs within the first encounter, whereas for 
b* < b < frscat it does not, but sufficient energy is radiated 
to put the binary into a bound state that eventually re- 
sults in a merger. A more precise definition arises in the 
context of the geodesic limit. The argument is that this 
threshold should gcnericaiiy be accompanied by behavior 
akin to zoom- whirl orbits in the geodesic limit [l9l. I20I |2~J] . 
For point particles orbiting BHs, the existence of zoom- 
whirl orbits is intimately related to that of unstable spher- 
ical orbits (at radii 3 < r/M < 6 for Schwarzschild BHs). 
One unusual property of these orbits is that an infinitesi- 
mal perturbation causes the particle to cither fall into the 
BH or "zoom" out on a bound elliptic orbit (if initially 
4 < r/M < 6 in the Schwarzschild example) or along a 
hyperbolic orbit (if initially 3 < r/M < 4). Conversely, 
a one-parameter family of geodesies that smoothly in- 
terpolates between capture by the BH (b < b*) and not 
(b > b*) will evolve arbitrarily close to one of these un- 
stable orbits as b — ► b* . The number of "whirls" close to 
the unstable orbit is given by 



C-rin|&-&* 



(1) 



where C is a family-dependent constant and T is inversely 
proportional to the Lyapunov instability exponent of the 
limiting spherical orbit. Although unstable spherical or- 
bits are a formal idealization, not realized in practice 
because of GW energy loss, BH mergers do indeed ap- 
proach a whirl-like configuration near the threshold of 
immediate merger. Such a configuration can in principle 
be sustained until all the excess kinetic energy, roughly 
equal to 2mi rr (7 — 1), is lost in GWs, where m,i rr is the 
irreducible mass of each black hole. This threshold is 



mine if it is consistent with the zoom-whirl picture, and 
if an expression of the form of Eq. |T]) holds. To this end 
we estimate the number of orbits n in the whirl phase 
in two ways: (a) using the puncture trajectories (n p ), 
and (b) using the GW flux measured far from the im- 
pact (new)- For method (a), in the scattering cases we 
define Tip cLS the total angle divided by 2tt traversed by 
the puncture from the initial position until it reaches a 
distance from the origin of twice the minimum distance 
(i.e., roughly twice the whirl radius). For merger cases, 
n p is counted in a like manner until the puncture crosses 
a distance 1/2 the whirl radius. Our current bracket 
gives b* /M = 3.35 ± 0.01, close to which we already see 
a "blurring" of the threshold, as binaries with b > b* 
do not separate to twice the whirl distance before merg- 
ing (see the inset of Fig. [T|). For method (b), we define 
now to be the number of GW cycles divided by 2 in the 
I = m = 2 component of the wave, from the initial time 
until the time when the I = 2 mode luminosity reaches 
1/2 its peak: see Fig. Q] for an example. We expect this 
estimate to be decent because the luminosity seems to 
be largest and roughly constant during the whirl phase. 
The "1/2 criterion" is somewhat arbitrary, but as long 
as it is applied consistently it should have little effect on 
our estimate of the slope T in Eq. |T]). 
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in parameter space. 



We now analyze sequence 1 in more detail to dcter- 



FIG. 1: An illustration of how we estimate hgw, here for 
b = 3.36M . The horizontal lines correspond to maximum, 
half-maximum and zero luminosity in the flux; the vertical 
line is when the flux drops below half maximum. For ref- 
erence, in the inset are puncture trajectories of one BH for 
this case (thick solid blue line), b = 3.34M (thin solid black 
line) and b = 3.39M (dash-dotted red line). All trajectories 
overlap early on during the whirl; the b = 3.39M trajectory 
peels off first before merging later, while the b = 3.34M case 
immediately spirals in after the whirl. The b — 3.36M case 
marginally peels off before merging, but it is clear from the 
luminosity that the merger/ringdown phase beginning after 
t/M ~ 280 is quite distinct from the preceeding whirl phase. 

The two estimates n p and tiqw are shown in Fig. [2] 
for immediate merger and nonprompt merger impact pa- 
rameters about b* . The plots indicate that a relationship 
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FIG. 2: Estimated number of orbits using puncture trajec- 
tories (left) and GWs (right) as a function of distance from 
b* for immediate and nonprompt merger cases in sequence 1. 
The uncertainties in the fitted slopes to Eq. (fTJ are purely 
from the linear regression analysis. 



FIG. 3: Total energy radiated (left) and final BH spin (right) 
vs. impact parameter from sequence 2, the latter calculated 
using several methods. The vertical dashed green (dotted 
red) line is the estimated immediate merger threshold b* (the 
scattering threshold & B cat). 



of the form ([I]) is valid, with a slope T w 0.2 to within 
~ 50% (ignoring systematic and computational errors we 
have not been able to account for) ; this is a factor of 2-3 
smaller than the analogous geodesic problem of a high- 
speed point particle scattering (prograde) off a Kerr BH 
with j 0.95. Given the large range of relevant T's in 
the geodesic case (cf. Fig. 9 of this provides reason- 
able evidence for zoom-whirl-like behavior in high-speed, 
comparable-mass collisions. 

The threshold of immediate merger b* appears to play 
a special role when we consider the amount of GW energy 
radiated and the final spin of the post-merger BH. To 
illustrate this point, in Fig. [3] we plot these quantities as 
functions of the impact parameter for sequence 2. 

For v ~ 0.75 the radiated energy increases by about 
one order of magnitude, from ~ 2.2% for b = to > 23% 
for b ~ b* . Two points are particularly noteworthy in 
this regard, (i) Even at this comparatively small boost 
v ~ 0.75 we comfortably exceed the maximum of 14±3 % 
reported for the ultrarelativistic limit of head-on colli- 
sions Grazing collisions with larger boosts, in turn, 
radiate enormous amounts of gravitational radiation: for 
one run of sequence 3 with v ~ 0.94 and b ~ b* the ra- 
diated energy is ~ 35 ± 5% of the CM energy, (ii) The 
maximum radiation as well as the maximum final spin 
(cf. below) is obtained near the threshold of immediate 
merger b* as opposed to b scat . The notion that excess 
kinetic energy drives zoom-whirl behavior seems to be 
consistent with the data, in that the maximum total en- 
ergy radiated near b* is approximately equal to the initial 
kinetic energy minus the spin energy of the final BH. 

These surprisingly large amounts of GW energy cor- 
respond to huge luminosities. Ref. showed that the 
high-energy, head-on collision of two nonspinning BHs 
could generate luminosities up to dE/dt ~ 0.01. For 
non head-on collisions and nonprompt mergers we ob- 
serve even higher luminosities. For instance, for v = 0.75 
and b ss b* the maximum luminosity is ~ 0.02, and 
extrapolation to v = 1 indicates that one might reach 
luminosities > 0.1, corresponding in physical units to 
~ 3.6 x 10 58 crgs _1 . This is the largest luminosity from a 
BH merger known to date, approaching in order of mag- 



TABLE I: Radiated energy, angular momentum and final spin 
for three representative sequence 2 runs. 

b/M E Tll d/M Jmd/M 2 Jrad/Jini jfin JQNM j AH 



0.00 


2.2% 


0.00 


0.00 


0.00 


0.00 


0.00 


2.74 


6.8% 


0.25 


0.25 


0.92 


0.91 


0.91 


3.04 


14.7% 


0.45 


0.40 


0.97 


0.96 


0.95 



nitude the universal limit dE/dt < 1 suggested by Dyson 

We conclude our analysis with a discussion of the fi- 
nal spin resulting from the merger. It has been argued in 
Washik et al. [23| that "no merger of equal-mass BHs can 
lead to a final BH with maximal spin parameter jg n « 1" , 
as long as the BHs are nonspinning. The maximum spin 
reported for the sequence studied in that work is 0.823. 
Larger final spins have been obtained from binary simu- 
lations with nonzero initial spin. Dain et al. [24j report 
the largest value measured so far (jg n = 0.922), although 
Fig. 3 in [25| might imply an even larger final spin. The 
latter work expresses doubts, however, whether merging 
equal-mass BHs can "produce a very-close-to-maximal fi- 
nal BH" . Our simulations suggest a different outcome for 
the merger of equal-mass, nonspinning BH binaries. In 
Table U we report the largest final spin measured so far 
in any numerical BH merger simulation. We further con- 
jecture that even equal-mass, nonspinning binaries can 
result in a final spin arbitrarily close to the Kerr limit 
j = 1. It is important in this context to bear in mind 
the difficulties in measuring the final spin with high ac- 
curacy (cf. [26]). These difficulties were our main reason 
to generate sequence 2. The high oscillation frequency of 
the ringdown signal as j — * 1 requires high resolution in 
the GW extraction zone. By using a smaller initial sep- 
aration, sequence 2 enables us to meet this requirement 
at tolerable computational cost. 

We have checked our results by calculating the final 
spin in a variety of ways: (i) we used energy balance ar- 
guments tO find jfin = Jfi n /Af| n = ( J ini - J rad )/ (M A DM ~ 

-Srad) 2 ; (ii) wc fit the quasinormal mode (QNM) fre- 
quency and damping time of the final BH and inverted 
them to obtain j'qnm (see e.g. [271]): (iii) we used the 
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equatorial circumfer ence of a Kerr BH C e = AirAI to find 
2ir Aah /Cg = 1 + y/l — wnere Aah is the AH area 
. These different estimates are compared in Table U 
for three selected impact parameters leading to merger, 
and shown in Fig. [3] for sequence 2 runs. Within our un- 
certainty estimates (~ 3 % for j'qnm, J ah and ~ 8 % for 
ifln) we observe good agreement throughout sequence 2. 
For 2.7 < b/M < b*/M we hnd j Rn > 0.9, and for b < b* 
our estimated final spins can be quite close to extrcmal- 
ity. For example, for b = 3.04M we directly measure 
j = 0.96 ± 0.03, and we expect further fine tuning of b to 
yield even larger values of j. Our estimates are substan- 
tially larger than the ones quoted by Shibata et al. fill ]. 
Given the difficulties in achieving the necessary numeri- 
cal accuracy, perhaps the apparent discrepancies are due 
to our increased resolution. 

IV. Conclusions. High-energy BH collisions are fertile 
ground for testing many ideas and conjectures in gen- 
eral relativity and high energy physics. We find that 
these collisions can radiate at least 35 ± 5% of the CM 



energy in GWs, that a merger can lead to the remnant 
BH spinning very close to extremal, and we display near- 
threshold phenomena akin to zoom-whirl in geodesies. 
We find no evidence of cosmic censorship violation. Our 
results are crucial for BH event generators in TeV-scale 
gravity: for four-dimensional spacetimes our results are 
consistent with Penrose's construction to within a fac- 
tor < 2. For the first time we show the dependence of 
the spin of the BH remnant on the impact parameter, a 
direct input in BH event generators. 
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